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SOME EFFECTS OF BLUNTNESS ON BOUNDARY-LAYER TRANSITION AND 
HEAT TRANSFER AT SUPERSONIC SPEEDS 1 

By W. E. Moeckel 


SUMMARY 

Large downstream movements oj transition observed when the 
leading edge oj a hollow cylinder or a jlat plate is slightly 
blunted are explained in terms oj the reduction in Reynolds 
number at the outer edge oj the boundary layer due to the de- 
tached shock wave. The magnitude oj this reduction is com- 
puted jor cones and wedges jor Mach numbers to 20. Concur- 
rent changes in outer-edge Mach number and temperature occur 
in the direction that would increase the stability oj the laminar 
boundary layer. 

The hypothesis is made that transition Reynolds number is 
substantially unchanged when a sharp leading edge or tip is 
blunted. This hypothesis leads to the conclusion that the 
downstream movement oj transition is inversely proportional to 
the ratio oj surface Reynolds number with blunted tip or lead- 
ing edge to surjace Reynolds number with sharp tip or leading 
edge. This conclusion is in good agreement with the hollow- 
cylinder result at Mach 3.1. 

Application oj this hypothesis to other Mach numbers yields 
the result that blunting the tip oj a slender cone or the leading 
edge oj a thin wedge should produce downstream movements oj 
transition by jactors ranging jrom 2 at Mach 3.0 to 30 at Mach 
15. The signijicance oj this result is discussed with regard to 
the possible reduction in over-all heat-transjer rate and jriction 
drag jor aircrajt jlying at high supersonic speeds. 

Mach number projiles near the surjaces oj blunted cones and 
wedges are computed jor an assumed shape oj the detached 
shock wave atjlight Mach numbers to 20. The dissipation and 
stability oj these projiles are discussed , and a method is de- 
scribed jor estimating the amount oj blunting reguired to produce 
the maximum possible downstream movement oj transition. 

INTRODUCTION 

In an investigation of the boundary layer on a hollow 
cylinder alined with the stream direction, it was discovered 
that the transition point moved downstream when the lead- 
I ing edge was slightly blunted (ref. 1). Similar results were 
obtained with a flat-plate whig in reference 2. A more ex- 
tensive investigation of the effects of leading-edge geometry 
on transition (ref. 3) confirmed previous results and led di- 
rectly to the explanation contained herein of the effect of 
blunting on transition. 

When a cone or wedge is blunted slightly (sketch (a)), the 
flow is changed in several ways, each of which could have a 

1 Supersedes NACA TN 3053, “Some Effects of Bluntness on Boundary-Layer Transition a 



noticeable effect on the transition location. A favorable 
static-pressure gradient is established near the vertex which 
could tend to stabilize the laminar layer. Downstream of 
the shoulder, however, the static-pressure gradient is ad- 
verse (for moderate supersonic speeds) because of the over- 
expansion around the shoulder and subsequent recompres- 
sion to the value corresponding to the unblunted cone or 
wedge. The effect of static-pressure gradient on transition 
is therefore inconclusive. 

In addition to the static-pressure gradient along the sur- 
face, the blunting produces a stagnation-pressure gradient 
normal to the surface. This gradient results from the var- 
iation in stagnation-pressure loss as the detached shock de- 
cays from the normal-shock strength at the vertex to the 
strength corresponding to the unblunted body at some dis- 
tance from the vertex. For inviscid flow, the stagnation 
pressure along each streamline remains constant downstream 
of the shock; hence, this gradient normal to the surface 
would persist along the entire length of the body. The 
stagnation-pressure gradient results in a shear layer whose 
thickness depends on the size of the blunted portion of the 
body. 

and Heat Transfer at Supersonic Speeds,’' by W. E. Moeckel, 1956. 


1 


2 


REPORT 1312 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


The fact that the entropy gradient produced by strongly 
curved shock waves might have appreciable effect on the 
development of the boundary layer is pointed out in refer- 
ences 4 and 5. Previously, the author of the present report 
had evaluated the shear profiles produced by detached shock 
waves near the surface of blunted flat plates. An explana- 
tion of the observed movement of transition in terms of 
these shear profiles was therefore sought. 

The interaction of the boundary layer with the shear 
profile produced by a detached shock wave is fundamentally 
a very difficult analytical problem, however, the condition 
of most interest is one for which the interaction of the two 
profiles is not important. Thus, if the shear profile produced 
by blunting is much thicker than the boundary layer, the 
rate of shear of the former is negligible compared with that 
of the latter. The boundary layer then develops in a region 
of negligible shear and in a layer whose Mach number is 
almost constant and less than that produced by a sharp 
cone or wedge. 

Of particular significance is the fact that the region of 
reduced Mach number near the surface is also a region of 
reduced Reynolds number. 2 Until the boundary layer en- 
gulfs this region, its stability and transition characteristics, 
as well as its friction and lieat-transfer characteristics, should 
be those associated with the reduced Reynolds number. 
This reduction in Reynolds number near the surface of a 
blunted body explains the downstream movement of transi- 
tion observed in references 1 to 3 and is the basis used in this 
report for comparing the boundary-layer characteristics of 
blunted and unblunted bodies. 

ANALYSIS 

The Mach number in the inviscid shear layer produced near 
the surface of a blunted cone or wedge increases continuously 
from the surface value to the value that would exist at the 
surface of the corresponding unblunted body. The Reynolds 
number per unit length at the outer edge of the boundary 
layer (outer-edge Reynolds number) therefore remains less 
than the free-stream (or unblunted) value until the boundary 
layer absorbs the entire shock-produced shear layer. If the 
transition point is determined primarily by the Reynolds 
number at the outer edge of the boundary layer, a progres- 
sive downstream movement of transition would therefore be 
expected as the leading edge or tip bluntness is gradually 
increased. The maximum downstream movement would be 
expected when the blunting is sufficiently great so that the 
outer-edge Reynolds number is close to the inviscid surface 
value for the entire laminar run. In the following sections, 
the maximum reduction in outer-edge Reynolds number is 
calculated, and the blunted area required to produce this 
maximum reduction over the entire laminar layer is estimated. 

REDUCTION IN SURFACE REYNOLDS NUMBER DUE TO BLUNTING 

At a station sufficiently far downstream of the vertex, 
where the surface static pressure for a blunted body is close 

2 This reduction in surface Reynolds number due to blunting and its effect on laminar heat- 
ing have recently been independently calculated in ref. 6 for hypersonic speeds. No at- 
tempt was made, however, to define the thickness and axial extent of the low Reynolds 
number layer or its effect on transition location. 


to that for the unblunted body, the Reynolds number near 
the surface can be written as 

Re n _ p n U n hi /p 

He i PiUi p n ytnMxPn 

where subscripts n and 1 refer to inviscid surface vaiues for 
the blunted and unblunted bodies, respectively. (All sym- 
bols are defined in appendix A.) These inviscid surface 
values will be assumed, as usual, to represent the outer-edge 
conditions that determine boundary-layer development. 

The use of Sutherland’s viscosity equation yields 


Re n ^t l \Vt n +S\M n 

R e \ \t n J Vi+*S/ Mi 


( 2 ) 


Dividing the numerator and denominator by the ambient 
static temperature t 0 and converting to Mach number 
functions yield 


where 


and 


H&n 712 -n -^n 

Re, 1 M, 


D= 


1+^Mf 
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2 
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2 


1 + y -~Ml to 
l+^M? to 


(3) 


(4) 


(5) 


The inviscid surface Mach number for the blunted body 
M n is determined by the ratio p, /P nj where p, is the static 
pressure at the surface of the unblunted body and P n is the 
stagnation pressure downstream of a normal shock at the 
flight Mach number M 0 . The inviscid surface Mach num- 
bers are shown in figure 1 as a function of M 0 for several cone 
and wedge angles. Since the total pressure P n is less than 
the total pressure at the surface for unblunted bodies, the 
surface Mach number M n for the blunted bodies is less than 
the surface Mach number for the unblunted bodies M The 
difference between M n and M, increases as flight Mach 
number increases. 

The Reynolds number ratio of equation (3) is plotted in 
figure 2 for the same cone and wedge angles as those in figure 
1. This ratio decreases rapidly as flight Mach number 
increases. If the transition Reynolds number is unchanged 
when the leading edge or tip is blunted, and if the blunting 
is adequate to cover the laminar boundary layer with a 
sufficiently thick layer of low Reynolds number air, then it 
should be possible to increase the length of laminar run by a 
factor inversely proportional to the Reynolds number ratio 
of figure 2. For slender cones and wedges, the possible in- 
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(a) Cone. 

(b) Wedge. 

Figure 1 . — In viscid surface Mach number for blunted cone and wedge. 


creases in laminar run range from factors of the order of 2.0 
at M 0 = 3.0, to 10 at M 0 = 8.0, and to 30 at M 0 =15.0. The 
significance of such large increases in laminar run for reducing 
the heat-transfer rate and friction drag for very high speed 
aircraft is apparent. 

Evidence that increases in laminar run of the magnitude 
indicated by figure 2 are actually attainable is presented in 
references 1 to 3. In reference 3, for example, the transition 
point at M 0 = 3.1 was moved downstream by a factor of 2 
(from 5 to 10 in. at a Reynolds number of 3.56X10 5 /in.) 
when the leading-edge thickness was increased from. 0.0008 
to 0.008 inch. This experimental movement of transition 
compares very favorably with the value 2.17 predicted on the 
basis of the Reynolds number reduction shown in figure 2(b). 

In reference 2, downstream, movements by factors ranging 
from 2.3 to 3.6 were observed for a blunted flat plate at 
various angles of attack at M 0 = 4.04. The movement pre- 
dicted by figure 2(b) for this Mach number is 3.57. For 
swept wings, little or no downstream, movement was observed 
in reference 2. This is in agreement with the expected weak- 
ening of the leading-edge shock due to sweepback. Whether 
downstream movements of the order of magnitude predicted 
by figure 2 are attainable at higher Mach number or for 
other body shapes remains to be established by further 
experiments. 



(a) Cone. 

(b) Wedge. 

Figure 2. — Effect of blunting on Reynolds number near surface of 
cone and wedge. Ambient static temperature, 392.4° R. 

estimation of bluntness required to obtain maximum 

MOVEMENT OF TRANSITION 

In order to determine the bluntness area required to cover 
the entire laminar boundary layer with a low Reynolds num- 
ber layer of negligible gradient, it is convenient to define a 
thickness of this layer which limits the Mach number to 
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values near the inviscid surface value. A suitable thickness 
is the distance from the surface to the streamline that passes 
through the sonic point of the detached shock wave (point 
where the Mach number just behind the shock is unity). 
From the vertex to the sonic point the stagnation pressure 
downstream of the shock does not vary greatly; consequently, 
the Mach number should remain near the inviscid surface 
value in the layer thus defined. 

An expression for the thickness will be derived under the 
assumption that the shear profile produced by the detached 
shock does not diffuse or dissipate, that is, the profile remains 
unchanged until it is engulfed by the boundary layer. The 
rate of dissipation of the shock-produced shear layer, which 
is discussed in appendix B, tends to increase the bluntness 
area required to produce a given thickness of the low Mach 
number layer. 

With dissipation neglected, the thickness of the low Mach 
number region can be estim ated by m eans of the detached - 
shock-wave theory of reference 7. In this theory, the de- 
tached shock wave is assumed to have a hyperbolic form in- 
dependent of the shape of the body that produces it. This 
form has been found to agree well with experimental results 
for a large range of body shapes in the moderate supersonic 
Mach number range (ref. 8) but becomes more questionable 
as Mach number increases. The location of the shock sonic 
point relative to the body is of particular significance for the 
present calculation. The assumption of reference 7 that the 
sonic line is straight and inclined normal to the mean flow 
direction is probably not a good approximation in the hyper- 
sonic speed range. Accordingly, the relative thickness of 
the low Reynolds number layer will be estimated from two 
separate assumptions: (1) the shock sonic point is located 
according to the method of reference 7, and (2) the ordinate 
of the shock sonic point is equal to half of the ordinate of the 
body sonic point. The latter assumption stems from con- 
sideration of a spherical-nosed body with the shock parallel 
to the surface at very high Mach numbers, where the dis- 
tance from, the shock to the body at the axis is of the order 
of 0.05 times the radius of the nose. The shock sonic point 
at very high Mach numbers occurs at a shock angle of ap- 
proximately 70° (corresponding to a radial angle of about 20° 
from the center of the spherical nose). The body sonic 
point, on the other hand, lies at a radial angle of about 45°. 
A layout of the parallel shock for these conditions yields 
approximately the result expressed by assumption (2). 

In order to estimate the thickness of the low Mach number 
layer, the Mach number in the layer is assumed to be con- 
stant at a value corresponding to the mass centroid of the 
layer. This Mach number, denoted by M c , is determined 
from the ratio pi/P c , where p { is the static pressure on the 
surface of the unblunted cone or wedge and P c is the total 
pressure downstream of the shock on the centroid stream- 
line. (A simple and satisfactory estimate of P c can be ob- 
tained by using the arithmetic mean of the stagnation pres- 
sures at the sonic point and at the vertex.) The continuity 
equation for the layer shown in sketch (b) can be written 

P c A n (A */A) Mc =PoA s (A*/A)„ o (6) 

where A n is the area of the low Mach number layer and A s 



is the free-stream. area of the stream tube between the axis 
and the shock sonic point S. If the bluntness of the body is 
defined as its cross-sectional area at the sonic point A SB 
(sketch (b)), the ratio of the area of the low Mach number 
layer to the blunted area becomes 

A n A s P 0 (A /A) M{) 

A sb A sb P c (A*/A) a , c 

For the assumption that the shock sonic point is located 
according to the method of reference 7 (assumption (1)), 
the ratio A S /A SB is given by 


p 

vvliere B=p^(A*/A) Mo and rj is the mean inclination of the 

sonic line defined in reference 7. For the assumption that 
the ordinate of the shock sonic point equals half the ordinate 
of the body sonic point (assumption (2)), A S /A SB is equal to 
0.5 for the two-dimensional case and 0.25 for the axisym- 
metric case. 

The area of the low Mach number layer defined by equa- 
tion (7) is shown in figure 3. For assumption (1), this area 
is seen to increase rapidly with increasing Mach numbers 
for the blunted flat plate (0^—0) and for the blunt-nosed 
cylinder (0 C =0). However, for wedge half-angles greater 
than 5° and cone half-angles greater than 10°, the area does 
not vary greatly with Mach number. For assumption (2), 
the area of the low Reynolds number layer is considerably 
smaller than for assumption (1). It is expected that the 
actual values will agree more closely with assumption (1) 
values at low supersonic speed and with assumption (2) 
values at very high speeds. Since the thickness of the low 
Reynolds number layer has been arbitrarily defined, it may 
be sufficiently accurate at this stage to state that the ratio 
A n /A SB is of the order of unity for slender blunted cones and 
wedges at all Mach numbers. 

With the thickness of the low Mach number and low 
Reynolds number layer thus defined, the blunting required 
to provide a low external stream Reynolds number for the 
entire laminar boundary layer to the expected or desired 
transition point can be estimated. This is done by calculat- 
ing the laminar boundary-layer thickness at the expected 
transition Reynolds number, which is based on conditions 
in the low Mach number layer near the surface. By equating 
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16 20 0 

Flight Mach number, Mq 
(a) Cone. (b) Wedge. 

Figure 3. — Area of low Mach number layer for blunted cone and 

wedge. 


this thickness to the thickness of the low Reynolds number 
layer, the required values of the bluntness area can be 
calculated. 

Thus, for blunted wedges the required ordinate at the 
body sonic point is 


Vsb 


(AJAsb) 


( 9 ) 


while for blunted cones (with 



27 * 1 , tr^tr 

( A n j A sb ) 



( 10 ) 


where AJA SB is given in figure 3 and r lt tT is the radius of the 
blunted cone at the expected transition point. Equations 
(9) and (10)- show that the amount of blunting required to 
cover the laminar boundary layer with a low Reynolds 
number layer is not large. For the wedge (eq. (9)), the ordi- 
nate of the body at the sonic point need be only of the order 
of magnitude of the boundary- layer thickness at the ex- 
pected transition point; for the cone (eq. (10)), the required 
radius of the body at the sonic point is of the order of the 
geometric mean of the body radius and the boundary-layer 
thickness at the expected transition point. 

The required bluntness is considerably reduced if the dis- 
placement effect is considered, since the low Reynolds number 
layer is moved away from the surface by an amount equal 
to the displacement thickness of the boundary layer (ref. 3). 
The required values of y SB should therefore be calculated 
with (8—8*) tr in place of 8 tr in equations (9) and (10). 
Expressions of 8 and 8 * for constant surface temperature 
were obtained from equations (18) and (22) of reference 9, 
based on the flat-plate theory of reference 10. The value of 
8 was assumed to correspond to a velocity ratio of 0.99. At 
the transition point, these expressions can be combined to 
yield (for 7= 1.40) 



- 3 - 48 a /£ 


(id 


where L tr is the distance along the surface to the transition 
point and C is the proportionality constant in the linear 
viscosity-temperature variation. For cones, this expression 
is divided by V3. In terms of (5— «*)*, the bluntness 
required to cover the laminar boundary layer with a low 
Reynolds number layer becomes 


for wedges and 


y SB ^ (8/L) tr (8*/L) tr 
Lit (AJAsb) 


( 12 ) 



mKHaj 


(A n /A SB ) 




(13) 


for cones. 

With equation (12), the calculated bluntness areas agree 
as closely as could be expected with the experimental values 
that produced the maximum downstream movement of 
transition in the experiments of reference 3. This maximum 
downstream movement was found to take place for a leading- 
edge thickness of about 0.008 inch, which is about two-thirds 
of the calculated value. Further increases in leading-edge 
thickness had no appreciable effect on transition location. 


INVISCID MACH NUMBER PROFILES FOR BLUNTED CONES AND WEDGES 

Although the maximum effect of blunting on boundary- 
layer development and transition depends on the portion 
of the shock-produced shear layer near the surface, the entire 
shear profile is of interest if the outer edge of the boundary 
layer moves out of the low Reynolds number layer defined 
in the preceding section. In order to determine the nature 
of the entire shock-produced shear profile, the shape and 
location of the shock must be prescribed. For moderate 
supersonic speeds, the hyperbolic form assumed in reference 
7 is adequate; but as the flight speed approaches the hyper- 
sonic range, the shape of the shock is increasingly influenced 
by body shape. This situation arises partly because the 
region between the shock and the body becomes smaller as 
M 0 increases; consequently, characteristics from portions 
of the body far downstream of the sonic point reach the 
shock before it has decayed to its asymptotic strength. In 
addition, the overexpansion near the shoulder of a slender 
blunted body, which takes place at lower speeds, gradually 
becomes an underexpansion at hypersonic speeds, that is, 
a Prandtl-Meyer expansion from the sonic point fails to 
reduce the pressure to, or below, the asymptotic static pres- 
sure. A rather long process of reflection of expansion waves 
between the shock and the body must, therefore, take place 
before the asymptotic pressure is reached on blunted cones 
or wedges. This consideration also affects the distance 
required to obtain the inviscid surface Mach numbers and 
Reynolds numbers calculated in the preceding sections. A 
more accurate evaluation of the effect of blunting would 
include the variation of outer-edge Mach number and Reyn- 
olds number along the entire body due to the pressure 
gradient. 
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For computing the shock-produced shear profile, the 
narrowing region between shock and body as flight speed 
increases introduces difficulties, in that no general shock 
shape is available beyond the sonic point and the asymptotic 
profile may be so far downstream as to have no practical 
significance. It was nevertheless felt to be worthwhile to 
compute these asymptotic profiles for a very high M 0 with 
the hyperbolic shock form of reference 7, if only for com- 
parison with more accurate future computations based on 
experimented shock forms or exact characteristic solutions 
for particular bodies. 

The computation method is presented in appendix C, and 
the resulting asymptotic inviscid shear profiles are shown in 
figure 4 for flight Mach numbers from 2 to 20. Indicated 
on each profile is the thickness of the low Mach number 
layer as defined in the preceding section. It is seen that 
this definition does, in fact, restrict the Mach number to 
values close to the surface value. 

The profiles for blunted wedges differ qualitatively from 
those of blunted cones at all Mach numbers. For the 
blunted wedges the Mach number gradient is zero at the 
surface; whereas, for blunted cones the gradient has a posi- 
tive value. This difference is traceable to the fact that, at 
the vertex of the shock, the flow area is proportional to the 
ordinate in the two-dimensional case and to the square of 
the ordinate in the axisymmetric case. 

The profiles of figure 4 should be good approximations for 
values of M 0 less than 5.0 but seem to become much too 



(a) Flight Mach number, 2.0; wedge, 
(b) Flight Mach number, 2.0; cone. 




(c) Flight Mach number, 3.0; wedge. 

(d) Flight Mach number, 3.0; cone. 



(e) Flight Mach number, 5.0; wedge, 
(f) Flight Mach number, 5.0; cone. 


Figure 4. — Inviscid Mach number profiles for blunted cone and wedge. 
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(g) Flight Mach number, 7.0; wedge, 
(h) Flight Mach number, 7.0; cone. 


5 * 





4 6 8 

Mach number, M 


(i) Flight Mach number, 10.0; wedge. 

(j) Flight Mach number, 10.0; cone. 



(k) Flight Mach number, 15.0; wedge. 
(1) Flight Mach number, 15.0; cone. 



(m) Flight Mach number, 20.0; wedge, 
(n) Flight Mach number, 20.0; cone. 


Figure 4. 


Concluded 
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thick at higher Mach numbers, particularly for the flat 
plate ( 0 ^= 0 ) and for the blunt-nosed cylinder ( 0 C =O). 
This thickness is associated with the overestimation of 
As/Asb by the method of reference 7 at these Mach numbers 
(see fig. 3). 

For higher cone and wedge angles, the shock decays more 
rapidly to the asymptotic strength, and the resulting pro- 
files appear to be more in harmony with expectations. 

Although the computed profiles are not reliable at high 
Mach numbers, they agree well with measured profiles at 
Mach 3.1 (ref. 3 ). If more accurate shock-produced profiles 
are desired for higher speeds, the shock form must be cal- 
culated for each body shape. Such computations would be 
useful for estimating the variation of transition location as 
the blunted area is gradually increased, but are not required 
for estimating the maximum downstream movement or the 
blunted area required to produce this movement. 

CHANGES IN LAMINAR RECOVERY TEMPERATURE, HEAT-TRANSFER 
RATE, AND FRICTION COEFFICIENT 

The downstream movement of transition due to blunting 
means that larger portions of the aircraft surfaces will be 
subjected to laminar, rather than turbulent, heat-transfer 
rates and friction coefficients. The blunting should, there- 
fore, produce substantia] reductions in over-all heat-transfer 
rate and friction drag. There is, however, an increase in 



Figure 5. — Effect of blunting on laminar heat-transfer rate, skin 
friction, and equilibrium temperature. 


laminar equilibrium recovery temperature corresponding to 
the reduction in Mach number and a change in laminar heat- 
transfer rate and friction coefficient due to the reduction in 
Reynolds number. These must be evaluated in order to 
estimate the magnitude of the advantages due to blunting. 

The heat-transfer coefficient and friction-drag equations 
of reference 11 are used for this estimate. Although these 
equations are based on the assumptions of constant specific 
heat and Prandtl number and no dissociation, they agree in 
order of magnitude with more exact numerical computations 
even at hypersonic speeds (ref. 12 ). 

The ratio of laminar heat-transfer rate with and without 
blunting can be written 


< 2 l ^1 \fe,l A 'Z>r 



(14) 


where D and r are defined by equations (4) arid (5), t e ,„ and 
t e i are equilibrium recovery temperatures with and without 

blunting, and g*(o, is the shear function of reference 

1 1 evaluated at the surface. This function is given in refer- 
ence 11 for several outer-edge Mach numbers M„ and 
several ratios ol surface temperature to outer-edge tem- 
perature <»/£„. 

The temperature-difference ratio in equation (14) can be 
written 
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t e, n t w ^0 

1+0.845 ( 

"y — V 

, 2 , 

l i 

to _ 

te, 1 t w t± 

to 

1+0.845 ( 

'7-1' 
, 2 , 

)m{ 

t U) 

to 


1 +0.2 M o | n , cn i\/T2\ 

T+oJrn {l+0 - mMn) ~T 0 


1 + 0.2 Ml 
1 + 0.2 M\ 


(1+0.169 M?)- 


1 1, 
to 


(15) 


where the laminar recovery factor is assumed to be 0.845 
and 7 is 1.4. 

The ratio of equilibrium surface temperatures with and 
without blunting is obtained from the definition of recovery 
factor: 

^^^V-r 1 ^ 0 - 845 (16) 

Jo — Jo n 


whence 


0.155 


1+0.2 M| 


+0.845 


0.155 


(17) 


1 + 0.2 M\ 


-0.845 


The laminar skin friction ratio is, from the equations of 
reference 11 , 

*{ c\ ^ w l\/f \ -M n h n 

(JM ,) » g _ V 0, + V Mi h 


(M„, 


T w, 1 


Z)V r 


g 


(0+M,) 


V5 


(18) 


The ratios of laminar recovery temperature, skin friction, 
and heat transfer for flat plates are shown in figure 5. 
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Although there is a slight increase in laminar equilibrium 
temperature for the blunted flat plate (this was observed 
experimentally in ref. 3), the laminar skin friction is reduced 
over the entire range of flight Mach numbers, and the heat- 
transfer rate is reduced except for wall temperature near 
equilibrium. (The rapid increase in the heat-transfer ratio 
near recovery temperature arises from the small increase in 
recovery temperature due to blunting. The heat transfer 
without blunting approaches zero for these values, whereas 
the heat transfer with blunting becomes small but is not yet 
zero.) Figure 5 shows that blunting the leading edge of a flat 
plate or cylinder can produce, in addition to the longer 
laminar run, a small but significant reduction in the skin 
friction and heat-transfer rate of the laminar boundary 
layer itself. 


COOLING REQUIREMENTS FOR STABILITY 

The static temperature at the edge of the boundary layer 
is considerably higher for a blunted cone or w^edge than for 
sharp bodies. The ratio tjti is, in fact, given by l/D (eq. 
(4)). This increase in outer-edge temperature means that, 
for a given surface temperature t w , the ratio tjt n is smaller 
than t w /ti. The outer-edge Mach number is also reduced. 
Shown in figure 6 are the outer-edge conditions for a blunted 
and unblunted flat plate and for a blunted and unblunted 
10° half angle cone for a surface-to-ambient temperature 
ratio of 4.0. These conditions are compared with two of the 
laminar stability limits given in references 12 and 13. This 
comparison shows that blunting moves the outer- edge con- 
ditions far into the stable region in the hypersonic speed 
range. (Although the stability-range curves shown are 
based on two-dimensional disturbance theory, recent com- 
putations by Dunn and Lin (ref. 13) indicate that three- 
dimensional disturbance theory also yields laminar stability 
to extremely high Reynolds numbers but that somewhat 
lower surface temperatures are required.) 

EFFECT OF BLUNTING ON HEAT-TRANSFER RATE NEAR THE NOSE 

In order to estimate more accurately the net decrease in 
heat-transfer rate due to blunting, it is necessary to deter- 
mine how the heat-transfer rate near the nose of the blunted 
cone differs from that on the pointed cone. An estimate of 
this difference can be made by comparing the heat-transfer 
rate for the sharp conical nose with that for the inscribed 
spherical nose (sketch (c)). 




Outer-edge Mach number, M m 

(a) 10° Half-angle cone. 

(b) Flat plate. 

Figure 6. — Effect of blunting on stability parameters for flat plate 
and 10° half-angle cone. Surface-to-ambient temperature ratio, 4.0. 


This heat- transfer ratio can be written as 


Qsp Asp A s p (Tq t A 

Q.c h c A c {Tg'l tv,) 


(19) 


where the subscripts sp and c refer to the spherical and 
conical noses, respectively. The area ratio of equation (19) is 


Asp = 9 tan Q c 
A c " cos e c 


(1— sin 6 C ) 


( 20 ) 


and the temperature difference ratio is 
T 0 t w 


TeA-tu 


y — 1 t 
^ *0 


1+^Ml 


( 21 ) 


— ii O+v^O-r 

1 +A_^M? V 2 7 to 


The mean heat-transfer coefficient for the spherical nose is 
assumed to be the stagnation point value presented in 
reference 14: 



( 22 ) 


where c=^-./— C p , st and (Nuf-yfRe^),,, is about 0.61 for 

a Prandtl number of 0.72 and for a ratio of wall temperature 
to stagnation temperature (tJT 0 ) close to zero (correspond- 
ing to cooled surfaces at a very high M 0 ). 3 The stagnation 
pressure coefficient C P , st is 1.84 for 7 = 1 . 4 . The mean cone 
heat-transfer coefficient is, from reference (15), 



(23) 


(c) 


3 Since the publication of ref. 14, Reshotko and Cohen have found that the expression for 
c given therein for supersonic flow is in error. The correct expression for this constant is that 
given above. 
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where He w ) c is 0.51. The ratio of heat-transfer co- 

efficients, therefore, becomes 



/ Nu \ 




p 

X 

xV H 6 w ) sp 

4 / Nu \ 

(: j£«°‘ »‘i 

L T 0 p c Y 

r'*" u pj 

(24) 


where P sl is the stagnation pressure behind a normal shock 
at Mach number M 0 . The heat-transfer ratio obtained by 
substituting equations (20), (21), and (24) into equation 
(19) is plotted in figure 7 for cone half-angles of 10° and 20° 
and for a surface-to-ambient static- temperature ratio of 
1.0. The over-all heat transfer for the inscribed spherical 
nose is seen to be less than that for the conical nose. The 
blunted nose, therefore, has the advantage of a lower heat- 
transfer rate near the vertex as well as along the downstream 
surfaces. 

DISCUSSION 

The preceding sections have shown that the Reynolds 
number per unit length at the outer edge of the boundary 
layer is lower for blunted fuselages and wings than for un- 
blunted ones. The limited data available agree with the 
conclusion that the transition location can be increased by a 
factor of the order of the ratio of the surface Reynolds 
number without blunting to the surface Reynolds number 
with blunting. This factor increases rapidly with increased 
flight speed, particularly for moderately slender wings and 
bodies. 

As an example of the magnitude of this effect, a 10° half- 
angle cone at a Mach number of 15 will be considered. If 
the transition point is located 1 foot downstream of the 
vertex without blunting, it might, on the basis of figure 2, 
be moved 25 feet downstream, of the vertex if the tip is 
blunted. 

The bluntness required is, from equations (11) and (13), 


^= 1.15 

Ptr 



where the value of A n IA SB corresponding to A S IA SB =0.25 

c 

has been used. If y — is of the order of 10~ 6 , then the re- 

quired value for ysB is about 11 inches. The ratio of the 
blunted area to the cross-sectional area of the cone at the 
transition point is, therefore, approximately 0.03. If the 
transition point (25 ft) is near the end of the body, the over- 
all heat-transfer rate would be reduced by blunting from 
the value corresponding to almost completely turbulent 
flow to the value corresponding to completely laminar flow. 
The blunted cone would, therefore, heat up much more 
slowly than the pointed cone and would require much less 
coolant to maintain a given surface temperature. The ratio 
L t rly n is about 600; therefore, the effect of dissipation of 
the shock-produced shear layer can probably be neglected 
(see appendix B). 

Furthermore, during the heating process the ratio of sur- 
face temperature to outer-edge temperature remains much 
lower for the blunted cone or wedge (fig. 6) so that the 
advantages of cooled surfaces with regard to laminar stabil- 
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Figure 7. — Comparison of laminar heat-transfer rates for spherical 
and conical noses. Surface-to-ambient static-temperature ratio, 1.0. 


ity prevail longer than for the pointed cone or sharp-edged 
wedge. Both the low surface Reynolds number and the 
higher outer-edge temperatures work toward preservation of 
laminar flow for a much larger distance along the surfaces of 
blunted bodies and wings. 

These advantages with regard to increased laminar run 
and increased laminar stability appear to involve no serious 
disadvantages. The friction drag is reduced, and the total 
drag should not increase appreciably for the small required 
values of the bluntness ratio. Reference 16 shows that, for 
spherical-tipped cones of fixed total length, the total drag to 
Mach number 7.0 is very near the value obtained for the 
sharp-tipped cone for ratios of nose diameter to maximum 
body diameter less than 0.25. 

The quantitative effects of blunting on transition location 
previously computed are based on the hypothesis that the 
transition Reynolds number is substantially unchanged when 
a body with a sharp tip is blunted. Although this hypoth- 
esis produces good agreement with the experimental results 
of reference 3, the possibility should certainly be kept in 
mind that at higher Mach numbers or with other body 
shapes the transition Reynolds number may be altered by 
such factors as pressure gradient and outer-edge Mach and 
Reynolds numbers. Furthermore, as the length of laminar 
run increases, the possibility of premature transition due to 
surface roughness or stream turbulence also increases, and 
the dissipation of the shock-produced shear profile becomes 
important. Whether any of these factors will seriously re- 
duce the attainable downstream movement of transition due 
to blunting remains to be determined experimentally. 

Many theoretical problems also require solution before 
the quantitative effects of blunting on transition can be 
predicted with confidence. One basic problem, of course, is 
that of the development of a laminar boundary layer in a 
nonuniform external stream. Solution of this problem 
would establish the magnitude of external shear that is 
negligible and, consequently, the conditions for which the 
boundary layer can be assumed to develop in a layer of 
reduced Reynolds number corresponding to the mean value 
near the surface. This solution might reveal whether, as 
indicated by the results of reference 3, it is sufficient, in 
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general, to obtain the maximum, reduction in Reynolds num- 
ber only for the inner half or two-thirds of the boundary 
layer rather than at the outer edge. The latter question, 
however, involves predicting the location of transition for 
various velocity profiles, which cannot as yet be done even 
for laminar layers in a uniform externa] stream. Since the 
required blunting is small, however, this question appears to 
be of secondary importance. 

CONCLUDING REMARKS 

It is clear from the preceding discussion that many ques- 
tions remain unanswered in this report. The principal ob- 
servation that the Reynolds number and Mach number near 
the surface are reduced by blunting and also the approxi- 
mate magnitude of the reductions are fairly well established. 
The assumption that the boundary-layer development should 
be determined primarily by the reduced Reynolds number 
and Mach number near the surface rather than by the flow 
outside the inviscid shear layer also seems reasonable. The 
principal benefits from blunting, however, lie in the hyper- , 


sonic speed range, where many of the quantitative results 
calculated herein are subject to corrections whose magnitudes 
are as 3 ^et unknown. Qualitative estimates indicate that 
some of these corrections, such as the displacement effect or 
the pressure gradients, either inviscid or self-induced by the 
boundary layer, should have a favorable effect on the down- 
stream movement of transition. Other effects, such as sur- 
face roughness, stream turbulence, or changes in transition 
Reynolds number, may tend to limit the downstream move- 
ment of transition to values less than those predicted. Dis- 
sociation at very high Mach numbers may have a significant 
effect on outer-edge conditions and, consequently, on the 
maximum transition movement to be expected. As usual, 
when so many unknown factors contribute to a phenomenon, 
experiment must be relied upon to determine which factors 
are dominant and which are of minor importance. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 21 , 1955 
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SYMBOLS 


area 

isentropic area contraction ratio from 
Mach number M to Mach number 1.0 
speed of sound 

^ (A*/A) Mq 

constant in linear viscosity-temperature 
relation 

stagnation-pressure coefficient 



1+0.2 M 2 n t x 
1+0.2 Ml t n 

shear function at surface (ref. 11) 

heat-transfer coefficient 
thermal conductivity of air 
length of conical tip 
Mach number 
Nusselt number 
stagnation pressure 
static pressure' 
heat-transfer rate 
Reynolds number 
recovery factor or radius 
cone radius at station where profiles are 
determined 

Sutherland's constant for air, 198.6° R 

stagnation temperature 

static temperature 

surface temperature 

velocity 

velocity downstream of normal shock 
ahead of spherical nose 
distance along surface 


y 

coordinate normal to surface 

P 


r 

t n ~i~S 

JL 


7 

ratio of specific heats, 1.40 

8 

boundary-layer thickness 

8 * 

boundary-layer displacement thickness 

V 

inclination of sonic line (ref. 7) 

6 C 

semivertex angle of cone 

6\v 

semivertex angle of wedge 

M 

coefficient of viscosity 

V 

m/p 

P 

density 

T 

shear force at surface 

<P 

shock angle 

Subscripts: 

C 

centroid 

c 

cone 

e 

equilibrium 

n 

inviscid surface values for blunted cones 


or wedges 

S 

sonic point on detached shock wave 

SB 

sonic point on body 

sp 

sphere 

St 

stagnation 

tr 

transition point 

W 

wedge 

w 

surface values 

0 

ambient conditions 

1 

inviscid surface values for unblunted 


cones or wedges 

00 

value at outer edge of boundary layer 


Superscript: 

' local conditions in inviscid shear layer 


12 


REPORT 1312 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


APPENDIX B 


DISSIPATION OF THE LOW MACH NUMBER LAYER 


In order to estimate more closely the amount of blunting 
required to maintain a given thickness of the low Mach 
number layer, the rate of dissipation of the shock-produced 
profile must be considered. The simplest method for esti- 
mating the rate is to consider the profile produced by the 
detached shock wave as a step function (sketch (d)) in which 



the outer velocity is that corresponding to the unblunted 
body and the inner velocity is that produced near the surface 
by blunting the vertex. The profile dissipation can then be 
considered identical to that at the interface of two parallel 
laminar jets emerging at the same static pressure. The 
equation for the velocity profile in the interaction region is 
given in reference 17 for the case when U x and U n differ by a 
small amount. The appropriate equation is 


U^Ui 

Un-Vi 




be obtained from the stability criterion for parallel jets 
developed in reference 18. This criterion states that the 


interface can become unstable if the quantity 


J) / du\ 
dyVby) 


vanishes in the interface profile. However, the profile is 
stable if this quantity vanishes only at points in the profile 
where the velocity satisfies one or both of the following 
inequalities: 

UKU.-a, 

U^>U n -{-a n 


where y n is the initial thickness of the low Mach number 
layer, Re Vn is Reynolds number based on y n and outer-flow 
conditions, and <£(<*) is the error function of a. Profiles 
calculated from equation (Bl) for Re 1/n =10 4 and for several 
values of x/y n are shown in figure 8. Apparently, the velocity 
near the surface does not change appreciably until x\y n is 
of the order of 1000. Although these profiles are valid only 
for small differences between Z7i and U n , the order of magni- 
tude of the dissipation remains the same for large differences 
(see fig. 4.11, ref. 17). The value of xVy n <1000 is, therefore, 
probably a good estimate for the length of run in which 
dissipation of the shock-produced shear profile can be 
neglected if this profile remains laminar. 

If transition to turbulence takes place in this layer, the 
length of run for which dissipation can be neglected is ap- 
preciably reduced. No experimental results are available 
to estimate under what conditions the shock-produced shear 
profile is likely to undergo transition. However, an indica- 
tion of whether transition is a possibility in this layer can 


These conditions assure that disturbances from either stream 
will not reach a layer in which amplification is possible. 

In terms of Mach number profiles, these conditions can be 
stated as follows: 

The profile is stable to two-dimensional disturbances if the 


quantity 


dM/dy 



-\ vanishes only at points where 


(B2) 


A/< My ~ 

V(7-l)(M 1 +2) 


or whore 


M> 7- 


1 


(y — 1)(2 — M n ) (B3) 

The latter condition cannot be satisfied for M n > 2.0. Con- 
sequently, the stability of the profile depends chiefly on 

whether — \ vanishes only where condition (B2) 


d y 


0 


7—1 


M 2 
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is fulfilled. Some sample computations based on the profiles 
of figure 4 indicate that condition (B2) is generally satisfied 
for the blunted-cone profiles but not for the blunted-wedge 
profiles. The latter profiles, therefore, are more inclined to 
undergo transition than the former. If transition occurs, the 
amount of bluntness required to produce a prescribed thick- 


ness of the low Mach number layer at a given station may 
be considerably greater than calculated on the basis of lami- 
nar flow. This discussion must necessarily be inconclusive, 
since the location or even the existence of transition cannot 
be established from stability theory alone. 


APPENDIX C 


MACH NUMBER PROFILES PRODUCED BY DETACHED SHOCK WAVES 


The Mach number profiles normal to the surfaces of cones 
and wedges with slightly blunted tips or leading edges can 
be calculated from the one-dimensional continuity equation 
if (1) the form and location of the detached shock wave are 
known, (2) the static pressure is constant normal to the 
surface, and (3) diffusion and dissipation of the profile are 
neglected. Condition (1) is most conveniently satisfied 
by using the detached-shock-wave theory of reference 7. 
Condition (2) is satisfied at stations sufficiently far down- 
stream of the nose or leading edge where the surface static 
pressure has reached, or closely approached, the value ob- 
tained with unblunted cones or wedges. At moderate 
supersonic Mach numbers, the required distance is of the 
order of 3 to 10 times the thickness of the blunted portion of 
the nose or leading edge. This condition is not quite satis- 
fied for blunted cones, because the flow' field approaches a 
conical distribution characterized by a gradual decrease of 
static pressure from the surface to the shock wave. But if 
the profile extends only a small portion of the distance from 
the surface to the shock w r ave, this gradient can be neglected 
without serious error. Condition (3) remains an assumption 
w T hose validity decreases as the distance along the body 
increases. It implies that the profile remains unchanged in 
form for an unlimited distance downstream of the vertex. 
As pointed out in appendix B, this assumption appears to 
be fairly good for distances of the order of 1000 times the 
thickness of the blunted portion of the body if the profile 
remains laminar. 

The profile computation is set up with the aid of sketch 
(e), which applies either for cones or wedges. If the sonic- 



point area of the body is used for reference purposes, the 
continuity equation can be written as 


or 


A' _ o / A \ 

A sb Jo P' (A»/4y d U s J 


(C2) 


where the primes refer to local conditions in the inviscid 
shear layer. Since the stagnation pressure along each 
streamline remains constant downstream of the shock, 
P'/Po is the stagnation-pressure ratio across the shock at 
the point where the streamline bounding the area A enters 
the shock. If the shock angle at this point is <p, then the 
total -pressure ratio can be written (ref. 19) 


^-1 

f 6Mg sin 2 ^ \» i 

< 6 V 

Po ' 

sin 2 <p+5/ ' 

v7 Ml sin 2 ip— 1/ 


(C3) 


The Mach number M' at the area A' can also be expressed 
in terms of shock angle by the relation: 


(C4) 


The function ( A*/A ) M < as well as P'/P 0 is a function of the 
shock angle <p at the point where the streamline crosses the 
shock wave. The differential d(^4/^4 sa ) of equation (C2) 
must be converted into a function of <p in order that the 
integration may be carried out from <p=90° to <p=<p h where 
<p i is the shock angle corresponding to the unblunted 
cone or wedge. The Mach number M' as a function of A' 
can then be obtained from equations (C4) and (C3). 

From reference 7 (eq. (5)), the relation between shock 
angle and shock ordinate for the assumed hyperbolic wave is 


where 


/ V V_ (xo/Vsb ) 2 
\VsB ) P 2 (0 2 tan 2 <p— 1) 


(C5) 


tan 2 **— 1 

VsB y S B 


(C6) 


and <p s is the shock angle at the shock sonic point. The ratio 
y sly sb is a function of M 0 and depends on whether the flow 
is tw T o-dimensional or axisymmetric. 

The area differential of equation (C2) can now be expressed 
as follows: 

For two-dimensional flow T : 


d ( — ^=d^ y \= b — r tan y sec2 y dy ~l 

\A S bJ \ysB ) p y.sB l_Q3 2 tan 2 
For axisymmetric flow: 


(C7) 


Af_ 

Asb 


PAiMBtfJ 

<A.\ 

Jo 

KAsbJ 


(Cl) 



, / Xp \ 2 r tan <p sec 2 <p dy ~l 

\y S B ) L(|8 2 tan 2 ip— l) 2 J 


(C8) 
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Combination of equations (C2) to (C8) yields the following final expressions for determining the variation of M' with A ' : 
For plane flow: 


A' 

Asb 


y' Vs iffTTZZi tPo C* ( 1 +0.2MI \» Mo r tan v secV 1, 


For axisymmetric flow: 


+=-2/3 2 (^Y(V tanVs— 1 

A S b \Vsb / \ 


\ Vo ( v / 1+0.2 Ml Y*M 0 1" tan <p sec fy "] i 

) Pi Jflo° \1 +0.2M' 2 / IT UP 2 tanV-l) 2 J ^ 


(C9) 


(CIO) 


In these equations, y f is the linear distance normal to the wedge or cone. For y'<^r u the area ratio A' / A SB in equation 
(10) is equal to 2 r x y' ly% B . 

Equations (C9) and (CIO) have been integrated numerically for several Mach numbers and for several wedge and cone 
angles. The resulting Mach number profiles are shown in figure 4. 
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